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Abstract 

We point out that in a large class of superstring-inspired E{Q) models, either an 
e"'"e~ collider operating at a center-of~mass energy ^/s = 1.5 TeV or higher must detect 
the pair production of charged or neutral exotic leptons, or an e'^e" collider with ^/s > 
300 GeV must discover at least one light neutral Higgs boson with invisible branching 
ratio exceeding 50%. If neither of these two signals is seen, the lightest neutral exotic 
lepton would overdose the universe, and the model could be completely excluded, 
independent of the values of the numerous free parameters. Future Higgs searches 
might lower the energy of the e'^e" collider needed to test these models decisively. 
The only assumptions we have to make are that -R— parity is exact, so that the lightest 
exotic lepton is stable if it is lighter than the lightest neutralino, and that no 50(10) 
singlet scalar gets a vacuum expectation of order lO^'' GeV or higher. If the second 
condition is violated, the model effectively reduces to an SO{ld) model as far as collider 
experiments are concerned. 



*Heisenberg Fellow 



1) Introduction 



Superstring-inspired £"(6) models [|l| contain a large number of particles in addition to 
those present in the Standard Model (SM): Superpartners of the known matter fermions 
and gauge bosons; scalar di- or lepto-quarks; extended gauge and Higgs sectors; and new 
"exotic" quarks and leptons. Indeed, this class of models can almost be considered to be a 
maximal (weakly interacting) extension of the SM. It is this aspect, rather than the by now 
quite tenuous connection to superstring theory 0, that keeps interest in these models alive 

i- 

Unfortunately most of the new particles predicted by E{6) models could be very heavy. 
In the absence of a comprehensive theory of supersymmetry breaking we are not able to give 
firm upper bounds on sparticle masses . The masses of the new gauge bosons contained in 
such models can be made very large by postulating large vacuum expectation values (vevs) 
for Higgs fields that are singlets under the SM gauge group. Finally, most of the exotic 
fermions reside in vector-like representations of the SM gauge group; they can therefore also 
be made very heavy. In the simplest of these models, where one requires gauge symmetry 
breaking to be triggered by radiative corrections P, |^, these possibly large scales are in 
fact all related: The vevs of most Higgs singlets cannot exceed the SUSY breaking scale 
significantly, and the same vevs also give rise to the masses of the exotic leptons, with 
Yukawa couplings of order 1 or less if the theory is to remain weakly interacting all the way 
up to the scale of grand unification (GUT). This allows to derive non-trivial relations 
between some of these masses, but unfortunately does not exclude the possibility that they 
are all very large. 

There are exceptions to this rule, however. Each fermion generation of -£"(6) contains 
27 degrees of freedom. Two of those are SM singlets, commonly called uji and A^. The 
resides in the 16 of 5*0(10), together with the 15 degrees of freedom that form a complete 
fermion generation in the SM. The fields might be exactly massless; alternatively, they 
might get very large masses through non-renormalizable operators if some scalar un field 
gets a vev at an "intermediate" scale around 10^° GeV or more, which could also give rise to 
see-saw type neutrino mass matrices [|^. Either way it is very difficult to derive significant 
constraints on these fields.[| 

In contrast, the superfields are singlets under 5*0(10). Their fermionic components 
can acquire masses at the weak scale only by mixing with the neutral components of exotic 
SU{2) doublet fermions |^; there are no terms of the type in the superpotential, since 

is not a singlet under the complete gauge group. Note that the vevs that give rise to 
this mixing also break the SU{2) x U{1)y symmetry of the SM, which means that they 
contribute to the masses of the W and Z bosons. This allows to derive a firm upper bound 
lOH of just over 100 GeV for the mass of the lightest eigenstate resulting from this mixing, if 



the relevant Yukawa couplings are required not to have a Landau pole below the GUT scale. 
This bound can only be avoided if some scalar has a vev of order 10^° GeV or more, and 
if the superfields containing the light exotics have non-renormalizable couplings to this vev. 
However, the same vev would also allow to give very large masses to all new gauge bosons 



^Successful nucleosynthesis in the early universe requires a large freeze-out temperature, i.e. large masses 
for gauge bosons coupling to lyji, if i^u is exactly massless 1^; however, direct experimental searches |^ by 
now constrain these gauge bosons to be heavier than several hundred GeV anyway. 
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and charged new matter fermions; at scales below (N) one would then just end up with 
the minimal supersymmetric standard model (MSSM) 0. We will therefore assume that 
scalars can only get vevs of the order of the weak or SUSY breaking scales, i.e. a few TeV 



or less; in this case the bound of ref.[jTO[ holds. 

The existence of an upper bound on the mass of an exotic fermion does not yet mean 
that we can test this model decisively, however. To begin with, present e^e~ colliders do not 
have sufficient energy to produce even the lightest exotics for all combinations of parameters. 
This will change as soon as the next generation of linear e~^e~ colliders goes into operation. 
However, even here the neutral exotics might be impossible to find. We will argue that 
the lightest neutral exotic lepton is likely to be stable. In this case it will only give rise 
to an observable signal at colliders if it is produced in association with a heavier exotic. 
Unfortunately no upper bound on the masses of these heavier states can as yet be given. 

In this seemingly hopeless situation cosmological arguments come to the rescue. It turns 
out that the upper bound on the mass of the lightest exotic decreases as the experimental 
lower bounds on the masses of the heavier exotics increase. Moreover, in this situation 
the light exotics are forced to have smaller and smaller SU{2) doublet components, which 
suppresses their couplings to gauge and Higgs bosons. Since the light exotics are stable, 
some fraction of such particles produced in the very early universe is still around today 
[ [Till . This fraction, and hence the contributions of such exotic Big Bang relics to the present 
matter density, is (approximately) inversely proportional to the annihilation cross section 
of the exotics. The crucial observation is that light, singlet -dominated exotics will have 
small annihilation cross sections by virtue of their suppressed couplings, and hence large 
relic densities. It is therefore clear that at some point the lower experimental bound on the 
masses of the heavier exotic leptons can force the cosmological relic density of the light exotic 
to become unacceptably large. We find that, before this happens, light neutral exotics will 
become accessible to Higgs boson decays, leading to a large invisible branching ratio for at 
least one light Higgs boson. It is the purpose of the present paper to study these connections 
quantitatively. Of course, the heavier exotics might well be discovered at rather low masses, 
in which case the relic density constraint could be fulfilled easily. In this optimistic scenario 
one can study the properties of the new fermions in order to further test the model. We wish 
to emphasize here that at some future time, the combination of collider searches and the relic 
density constraint will test the model decisively, independent of the values of the (many) 
free parameters. Unfortunately we find that a 500 GeV e^e~ collider is not quite sufficient 
to probe the entire parameter space; one will have to push the energy to 1 or even 1.5 TeV 
in order to close the last loophole. Nevertheless we find it remarkable that a decisive test is 
possible at all; after all, most experimental bounds can be evaded by simply increasing the 
"new physics" scale. 

The remainder of this paper is organized as follows. In Sec. 2 we describe those parts of 
the model that are of relevance to us, i.e. the exotic leptons and Higgs bosons; we also give 
the couplings of the (neutral) exotics to gauge and Higgs bosons. In Sec. 3 the cross sections 
for the production of neutral exotics at e~^e~ colliders are listed. In Sec. 4 the calculation 
of the cosmological relic density of the light neutral exotics is discussed; some care must be 
taken here, due to the prominent role played by s— channel exchange diagrams, as well as 
the presence of a second light exotic, which can suppress the relic density by co-annihilation. 
Sec. 5 contains our numerical results, and Sec. 6 is devoted to a brief summary. 
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2) The Model 



For defineteness we will work in the framework of the minimal rank 5 subgroup of -^(6) [|T|, 
which does not require the introduction of an intermediate scale between the GUT and weak 
scales. However, our analysis would go through with only minor modifications in models 
with a gauge group of rank 6 at the weak scale, as well as in models where some field has 
a vev around 10^" GeV or more. As already stated in the Introduction, the only important 
assumption we have to make is that none of the scalar N fields gets a vev much above 10^ 
GeV. 

In E{6) models each 27-dimensional representation contains one superfield, which 
is an SM singlet, as well as the exotic SU{2) doublet superfields H and H, whose scalar 
components have just the right quantum numbers to serve as the Higgs bosons of the MSSM, 
i.e. to provide masses for the W and Z bosons. At least one of the A^ scalars also has to 
get a vev in order to give a sufficiently large mass to the single new Z' boson of the rank-5 
model. Using rotations in generation space, we can always work in a basis where only one 
H^, one and one A^ field have non- vanishing vev; following the notation of ref.0 we 
call these true Higgs fields H^, and A^s. Their fermionic superpartners then mix with 
the superpartners of the three neutral gauge bosons to form the six neutralino states of this 
model §. 

Here we are interested in the first two generations of fermionic H, H and A^ fields. They 
obtain their masses from the superpotential 



lep 



i,j,k=l 



HiHjNk. 



(1) 



In order to give masses to all charged exotic lepton^, we at least have to allow those couplings 
where exactly one of the indices i,j,k in eq. ([^) equals 3, the other two being either 1 or 2. By 
further rotations between fields of the first two generations we can define A123 = A213 = 0, 
without loss of generality; the contribution of these couplings to the mass matrix of the 
charged exotic fermions is then diagonal, with 



mj± 



^223^, 



(2a) 
(2b) 



where we have introduced x = (N^)- 

In the basis [H^, H^, Ni, H2, H2, N2), the contribution of these couplings to the mass 
matrix of the neutral exotic leptons reads: 



M 



LO 







Al31^^ 








Al32^^ \ 







Asiit^ 








A312^^ 


Al31^^ 


A311^^ 





A231^^ 


A321^^ 











A231^^ 





^2 


A232^^ 








A321W 







A322^^ 


V Ai32t) 


A312^^ 





>^232V 


A322^^ 


J 



(3) 



*We call these fields "leptons" merely in order to indicate that they do not have strong interactions; they 
need not carry the same lepton number as the charged leptons and neutrinos of the SM 
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with V = (H^), V = (-f^a). The neutral exotic leptons are Majorana fermions; the mass 
matrix is therefore symmetric. 

The couphng A333 in eq.(|l|) must also be non-zero, in order to avoid the presence of a 
dangerous axion (see below). Terms where none of the three subscripts equals 3 do not 
contribute to the mass matrices. Finally, there could be contributions to the superpotential 
(|I]) where only one of the three indices is not equal to 3. In this case the charged exotic 
leptons would mix with the charginos, and the neutral exotic leptons would mix with the 
neutralinos. However, if such couplings are present, the rotatations in field space that define 
the basis where only H^, and have non-zero vevs would depend on the renormalization 
scale. Worse, all Hi and Hi would then couple to SM quarks and leptons at least at the 
one loop level, which could lead to dangerous flavour changing neutral currents (FCNC) 
1^ . We therefore forbid these terms. Fortunately, this is not only technically natural, but 



even follows automatically if one requires all potentially dangerous terms in the low energy 
superpotential (which can lead to fast proton decay, large neutrino masses, or tree-level 
FCNC) to be forbidden by a single discrete Z2 symmetry |T^ . 



This further restriction of the allowed terms in eq.([l|) greatly simplifies the calculation, 
and leads to a much more appealing model, but has no significant impact on the properties 
that are of interest to us. This can be seen from the fact that even if all terms in eq.(0) 
were present, in the limit — > the 12 x 12 neutral fermion mass matrix would have 
two zero eigenvalues, corresponding to SM singlet fermions A^i, N2- Since all couplings in 
eq.([^) are required to be of order 1 or less H, [l^, this observation implies that even in the 



realistic situation with non-vanishing v and there will be two neutral Majorana fermions 
whose masses are roughly of order Mz or less. Moreover, the upper bound on the masses 
of these neutral fermions will decrease as the (experimental lower bounds on the) masses of 
the charged exotic fermions is increased. We repeat, these crucial properties of the model 
follow from eq.(|I]) without further assumptions. For the remainder of our analysis we will 
concentrate on the simpler case where the mass matrix (0) is decoupled from the neutralino 
mass matrix, just to avoid needless complications. 

A little calculation shows that the determinant of the matrix (0) is proportional to {vvxY\ 
this means that all three vevs have to be non-zero if all eigenvalues are to be non-vanishing. 
It is also quite easy to see that in the limit '^l*^ ^ v,v, and with all couplings 0{1) or 
less, there will be four large eigenvalues, approximately equal to ±^«^± and ±m]^±. The 
other two eigenvalues are then of order A^ft;/m^± . This is why the upper bound on the 
smaller eigenvalues decreases with increasing mass of the charged exotics, as stated above. 



Indeed, in ref. it was shown that the mass m^o of the lightest neutral exotic lepton is 
maximized if all entries of the mass matrix (|^) are of the same order (or zero). Numerically 
one has m^o < 110 GeV if 'm^± > 45 GeV. 

In general the matrix (|^) has to be diagonalized numerically. The resulting eigenstates 
L° are given by 

4° = Ef/..iv„ (4) 

where we have introduced the 6-component vector = {Hi, Hi, Ni, H2, H2, N2), and U is 
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an othogonalC] matrix chosen such that 

UMloU^ = diag(m^o) 



1,...,6. 



(5) 



The couphngs of the charged and neutral exotic mass eigenstates to the standard Z boson 
are then given by the Lagrangean 



ZLL 



2cos6'i4/ 



^I^7'^L±(1 - 2sin2^H^) + ^ 4V75^,% 

a=l *ii=l 



(6) 



where g is the SU (2) gauge couphng, 6'vi/ the weak mixing angle, and 



UaU. 



(7) 



A few comments are in order here. First, when writing eqs.(^),@ we have allowed 
the masses of the neutral exotic leptons to have either sign. We could also insist that all 
m^o > by inserting appropriate factors of i in the matrix U, which would then no longer 

i 

be orthogonal (but would still be unitary, of course). In this case the diagonal couplings 
of the Z boson to neutral exotics would still be purely axial vector, but the off-diagonal 
couplings would become purely vector if the two corresponding eigenvalues have opposite 
signs. Secondly, recall that the L° are (4-component) Majorana spinors. This means that 
their couplings appearing in Feynman diagrams are twice as large as those in the Lagrangean 
(^. Finally, following ref.|^, and contrary to the usual MSSM notation, we have defined 
/7i.2,3 to have hypercharge +1/2, while the hypercharge of -f^i,2,3 is —1/2; this explains the 
overall sign of the coupling tij in eq.(0). 

In our calculation of cosmological relic densities we also have to specify the Higgs sector 
of the model. Recall that we are working in a basis where only H^, and have non-zero 



vevs. The relevant part of the Higgs potential is then given by |15 



Higgs 



r0|2 



+ A 



333 



(8) 



+ 



9^ 



8 cos^ Ow 



I //Op 



where g' is the U{1) gauge coupling. We will assume the Higgs mass parameters as well as 
the soft breaking parameter ^333 to be real. The vevs can then all chosen to be real |]16| , 
and the mass matrices for the real (scalar) and imaginary (pseudoscalar) parts of the neutral 
Higgs fields decouple. The former is given by P 



M 



(2\ 

V (2-^333 



25 ,2 



V 



A 



vx 



2 ^ 18- _ 
333 - 4 - + 

- ^g'A VX + Av 



(2A2 - 5 .'2^ 



2 

333 



2\l,,-'ig'^]vx + Av\ 

^^333 



-^g^]vx + Av 



j^g ] vx + Av 
Wx- - If 
(9) 



^We assume that all couplings in eq.(|l|) are real. 
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with A = X333A333. In writing eq.(|^) we have used the requirement that the first derivatives 
of the potential with respect to the fields should vanish in the minimum. Notice that 
the smallest eigenvalue m^o of the matrix (P) is again only of order f ^ + v'^, not of order of 
the Z' boson mass or the SUSY breaking scale |[T7[|: 



mlo < Ml 



cos'm + '^^l^^'^'f^^ sin2(2/5) + (4 sin^ /3 + cos^ /5 

9 9 ^ 



(10) 



where tan/3 = v/v (> 1). Eqs.(^ and (|TOD only hold at tree level. There are sizable radiative 
corrections from top and stop loops, as well as possibly from loops involving exotic (s) quarks 
and (s)leptons |T^. However, for our purpose their effect can largely be mimicked by allowing 
A333 at the weak scale to be quite large. The reason is that these corrections tend to increase 
the mass of the lightest scalar Higgs boson; eq. (0) shows that increasing A333 has the same 
effect. Moreover, the radiative corrections are only important if they involve large Yukawa 
couplings, e.g. that of the top quark. However, the introduction of these large couplings 
also reduces the upper bound on A333 which follows from the requirement that the model 
remains weakly interacting up to the GUT scale. The final upper bound on uij^o is therefore 
not much changed by these corrections. 

Since the model contains two massive neutral gauge bosons, only one physical pseu- 
doscalar A survives. It can be written as p 



1 1 ^ f ^-n VV 

V 



A = A.^ { -m + -HI + ^iV3 



Jv/v + v/v + vv/x^ 

V2^ yv 



^^iiPj^iH,), (11) 

where we have introduced the complex 3-component vector H = (^Hl,Hl,N3j, and the 
symbol 3 denotes the imaginary part. The mass of the physical pseudoscalar is given by 

o . , V VV\ , , 

ml = -\3s3A333X - + - + ^ • (12) 

\V V X'^ J 

Notice that this state would be massless if A333 = 0, as mentioned earlier. There is also a 
physical charged Higgs boson, with mass 

<± = M^(l - ^) + mi + A333A333-. (13) 

Note that this Higgs boson can be lighter than the W boson if m^ is small and A333 is sizable. 

In general the mass matrix (y) is most easily diagonalized numerically. Since it is real 
and symmetric, the diagonalization can be achieved by an orthogonal matrix 5*: 

SM\S^ = diag{mlo,mlo,mlo). (14) 
The physical scalars are then given by 

^° = 4^ E 5^.^ (Hj) . (15) 
V 2 j=i 
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The matrix S and the vector P introduced in eq.(|TT]) determine the couphngs of the physical 
Higgs bosons to SM particles. We need the Z — Z — h^ couplings in order to interpret bounds 
on the Higgs sector from searches at LEP; they are given by 



9zzh° 



gMz 



cos 9w 

Bounds on Z h^A decays involve the couplings 
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(S'iisin/? + Si2Cos(3) . 



9zAhO 



2 cose 



{SiiPi — 5*42 -P2 J 



(16) 



(17) 



w 



In Sec. 4 we will also need the couplings of the Higgs bosons to the quarks and leptons of 
the SM. We write the corresponding Lagrangean as 



C 



Hff 



9 



cos9 



w 



-f 



.if) I 



-^75)a + E4^^^° 



/, 



with 



.(/) 



d 



if) 



nif J 


f Pi/sin/5, 


f = 


u 




t 


2Mz 1 


[ P2/COS/3, 


f = 


e 




T, d, s, b 


nif 1 


Sii/sinP, 


f = 


M, 


c, 


t 


2Mz 1 


Si2/cos(3, 


f = 


e, 




T, d, s, b 



(18) 

(19a) 
(19b) 



Finally, we will need the couplings of the physical Higgs bosons to the neutral exotic 
leptons. They can most easily be expressed in terms of the couplings A^^-^, between the lepton 
and Higgs current eigenstates: 



'NNH 



j:K^,m{i-^,)N,H', + h.c., 

i,j=l k=l 



(20) 



where is the 6-component vector of Majorana spinors introduced in eq.(^, and H the 
3-component vector of complex neutral Higgs fields defined in eq. ([TT|) . The couplings X'^jj. 
are determined by the superpotential (|l]): 



A' 



123 



A' 



■^113) -^132 — -^131) -^153 ~ ^1237 -^162 — -^132, A231 — A311, A243 — A213, 
261 ~ -^312, A342 = A23I, X'^5i = A321, A453 = A223! -^462 = -^232, X'^Ql = A322; (21) 



this has to be symmetrized in the first two indices (A-^^ = X'jik), and all other A' couplings 
vanish. Recall that without loss of generality we can work in a basis where the mass matrix 



for the charged exotic leptons is diagonal, i.e. A123 = A213 = 0; this implies \\ 



153 



A' 



243 





in this basis. The interaction between lepton and Higgs mass states can then be written as 

1 6 3 3 6 



^V^iJ=lk=l 1=1 m,n=l 

A E E K.kPkA E f/™^nA75l^°, 

^V^ij=lfc=l m,n=l 



+ 



(22) 
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where the orthogonal matrices U and S have been defined in eqs.(|^) and ([T^), respectively, 
and the eigenvector P of the pseudoscalar mass matrix is given in eq.(0). Recall that the 
are Majorana fermions, and the physical Higgs bosons and A are described by real 
fields. 



3) Constraints from e^e Colliders 

In this section we discuss how new physics searches at existing and future e^e^ colliders 
can constrain the class of E[Q) models we are considering. At present the most stringent 
constraints on both the exotic leptons and the Higgs sector come from LEP. To begin with, 
the failure to observe the charged exotic leptons immediately implies that their masses must 
exceed 45 GeV, since eq.(^ shows that they always couple with essentially full gauge strength 
to the Z boson. 

The decay of the Z into two neutral exotic leptons needs a somewhat more detailed 
discussion. The heavier L° states can always decay into plus a real or virtual Z or 
Higgs boson. However, the lightest state can only decay into SM particles or the right- 
handed neutrino state un, or their superpartners. Such decays could be due to terms in the 
superpotential of the type 111,21 l^r, -ffi,2Q'LC?i? or Hi,2(1lUr, where qi and II are quark and 
lepton SU{2) doublet superfields and cr, cIr and ur the corresponding singlets. However, 
couplings of this kind will in general generate tree-level FCNC; these terms are therefore 
severely constrained |T^, which is why they are often completely forbidden in (potentially) 
realistic models, e.g. by means of a discrete symmetry The model in general also allows 
couplings of the kind Hi 21l^r. Such couplings can lead to fiavor changing processes in the 
lepton sector (/x — >• ej, r /i7, etc), but only at one-loop level; they are therefore somewhat 



less tightly constrained |T9[ . If any of these couplings exists, the exotic leptons become odd 
under R parity. This means that the decay product of must contain a superparticle. Since 
L° is quite light, we can expect that only final states containing the lightest sparticle, which 
in almost all cases is the lightest neutralino Zi, can occur. The least tightly constrained of 
the above couplings, Hi 21li^r, will then lead to an invisible final state, if we assume that 
R parity is conserved. Of course, it is not at all clear that such decays are kinematically 
allowed, since Zi could very well be heavier than Indeed, in view of the fact that present 
bounds p, 1^ require the Z' boson mass to exceed several hundred GeV, we need a rather 
high value of the SUSY breaking scale in this model |^. We therefore conclude that L° 
most probably is either absolutely stable, or decays invisibly; in particular, this is always 
true if tree-level FCNC are suppressed by a simple symmetry and R parity is conserved. In 
this case L° cannot be detected directly by collider experiments. Since we want to devise 
a test of the model that works even under the least favorable circumstances we will assume 
that L5 is indeed invisible. 



The partial width oi Z ^ ^i-^^ decays is given by p!4| : 
TiZ ^ L°L°) = J^(2 - 5,,) ( 
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2M7 — m,- — m,- — bmirrin — ^ — ^ 



Ml 



2- 

' (23) 



where we have used the shorthand notation = m^^o. The couphng tij has been defined in 
eq.(0), and the 3-momentum in the Z rest frame |fc| is given by 



1^1 = ^ V (m| - mf - m|) ' - 4m2m2 . (24) 

Recall that in our convention the neutral lepton masses, and hence the bi-linear term in 
eq.(|23|), can have either sign. 

We have argued above that the combination i = j = 1 only contributes to the invisible 
width of the Z boson. However, even though the heavier neutral leptons will almost always 
decay inside the detector, their production might not lead to an experimentally observable 
final state if the mass difference to the lightest lepton is too small. We have quite conser- 
vatively required |mj| — \mi\ > 3 GeV for an "experimentally visible" L?. The signature 
for the production of these visible states is quite similar to that for the production of the 
heavier neutralino states of the MSSM. Searches for such neutralinos have so far not 



been sucessful; the bound on the resulting branching ratio is of the order of (a few times) 
10~^, the precise value depending on the masses of the particles involved. We have therefore 
imposed the upper bound 

r{Z visible exotics) < 0.13 MeV, (25) 

which corresponds to a bound on the branching ratio of about 5 ■ 10~^. 

The best bound on non-SM contributions to invisible Z decays comes from the measure- 
ment of rinv(Z)/r(Z this can be interpreted pl| as a measurement of the number 



of SM neutrino species, Ny = 2.988 ± 0.023. Taking into account that the model predicts 
Ni, > 3, this translates into the bound 



AFinv < 6.7 MeV (26) 

at 95% confidence level. Note that "invisible" includes everything that is not counted in any 
explicitly reconstructed final state; in particular, L^L^ final states where neither of the two 
leptons is at least 3 GeV heavier than L° are included here. 

At future high energy e~^e~ colliders the neutral exotic leptons can be produced via the 
exchange of a virtual Z boson. The cross section is given by 

a(e-e-^LM = ^(2-V^ ' ^^^"^"^^ 



E,E^-niinij + ^\k\'^^, (27) 

with Og = —1/4 and fg = 1/4 — siv? 6w \k\ is given by eq.(plD with Mz replaced by 

the center-of-mass energy y^, and E^ = \l mf + In principle this cross section also 
receives contributions from Z' boson exchange. However, if these contributions are sizable. 
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Z' exchange will also lead to observable effects in the pair production of SM fermions. Future 
e^e~ colliders will therefore be able to find evidence for a Z' boson with mass up to several 
times a/s. The discovery of a Z' signal would not only (obviously) rule out the SM, but would 
also allow significant new tests of the type of model we are studying here. In order to be 
conservative we therefore always assume that the contributions from Z' exchange are too 
small to be detectable. Finally, we will always assume that future e^e~ colliders will probe 
for the existence of charged exotic leptons almost up to the kinematical limit, since they 
have full strength couplings to both the Z boson and the photon. 

The Higgs sector of the model is also constrained by unsuccessful new physics searches 



at LEP. We have used the results of the ALEPH collaboration |^, which has published 
bounds on the ZZh^ and ZAh^ couplings of light Higgs bosons. Their numerical bounds 
can be approximated by: 



{Sii sin/5 + Si2 cos/3)^ < < 



0.025, my^o < 10 GeV 

O.OOSm^o - 0.025, 10 GeV < m;,o < 20 GeV 

0.0175m';,o - 0.275, 20 GeV < m^^o < 30 GeV 

0.025m^o - 0.5, 30 GeV < m^o < 60 GeV 



(28a) 



^ . 0.1, m/jO + niA < 81 GeV 

{SaPi - <{ Q / + - 81 GeV) + 0.1, m^o + < ^^^^^ 



where all masses are in GeV. 

As emphasized in Sec. 2, at least one of the three neutral scalar Higgs bosons of the model 
must have mass below 150 GeV or so; moreover, this boson will couple to the Z with full 
strength if the bound (p!0| ) on its mass is saturated. For some combinations of parameters 
the lightest Higgs scalar will have a very weak coupling to the Z; however, in this case the 
next-to-lightest Higgs scalar will have unsuppressed coupling to the Z, and its mass will 
also satisfy the bound (|l3)- In fact, the scalar Higgs sector of our model is quite similar to 



that of the MSSM with additional Higgs singlet superfield, where it has been shown |2^ that 
an e"'"e~ collider with -^i > 300 GeV has to detect at least one neutral scalar Higgs boson. 
Unfortunately this may not be sufficient to allow a significant test of the model, since this 
Higgs boson might look very similar to the single Higgs boson of the SM. Even if several 
Higgs bosons are discovered, it might be quite difficult to distinguish between the £"(6) model 
were are discussing here and the MSSM, since large x means that one of the three neutral 
Higgs scalars is quite heavy and essentially a pure singlet; the singlet components of the other 
two scalar Higgs bosons and the single pseudoscalar Higgs boson are then very small. We 
therefore try to avoid making assumptions about searches for neutral Higgs bosons at future 
colliders as much as possible. There is one exception, however: The light Higgs boson(s) of 
the model might have large branching ratios into neutral exotic leptons, which could lead to 
a large invisible branching ratio for these Higgs particles. This would be fairly distinctive; 
in the MSSM a large invisible branching ratio of a Higgs boson would imply that the light 
chargino will be discovered at LEP2 |2^. In contrast, we will always require that the charged 



exotic leptons are too heavy to be produced at the e+e collider under consideration. 
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The partial widths for Higgs decays into exotic lepton can be written as 



m'jj — {rrii ± mj)^ 



(29) 

As usual, the lepton masses rriij can have either sign here. The "+" ("— ") sign in the last 
term of eq.(]29|) applies if H is scalar (pseudoscalar), and A is the relevant HL^L^ coupling 
in the Lagrangean (p2D; recall that a factor of 2A appears in the relevant Feynman rule, due 
to the Majorana nature of When computing the total decay widths of the neutral Higgs 
bosons of our model, we include decays into SM fermion pairs as well as (for the scalar Higgs 
bosons) W^W~ and ZZ final states, where we allow one (but not both) of the gauge bosons 
to be off-shell. Expressions for the relevant partial widths can, e.g., be found in the recent 
review article [Q; of course, we have to use the couplings listed in eqs.(p^-(p^ here. Note 
that the quark masses in eqs.(|l^) are meant to be running masses, taken at scale niH] we 
also include the leading non-logarithmic QCD corrections |^ (which remain finite in the 
limit niq/mH — > 0) to the cc and bb partial widths. An accurate estimate of the decay widths 
of the light Higgs bosons of the model is necessary not only for the evaluation of branching 
ratios for exotic Higgs decays, but also for the calculation of the cosmological relic density, 
to which we turn next. 



4) The Cosmological Relic Density 

We argued in the previous section that the exotic leptons will have odd R parity if they 
have any Yukawa interactions with ordinary matter. The lightest exotic lepton will then 
be absolutely stable if it is lighter than the lightest neutralino Zi. The only way out would 
be to introduce R parity breaking interactions in the superpotential. However, in this case 
the model, although very complicated, would no longer be able to explain the existence of 
cold dark matter (CDM); this seems to be required for a successful fit of all data on 
large scale structure in the universe. The fact that Zi makes an excellent CDM candidate is 
often viewed as one of the strengths of the MSSM; a significantly more complicated model 
that gives up this advantage does not look very appealing to us. We will therefore require 
R parity to be conserved. 



Unfortunately the usual approximation for the calculation of the relic density is often 



not reliable in this model. Two of the three cases discussed in ref . as examples where the 
usual method is not applicable can occur here, possibly simultaneously: Narrow s— channel 
(Higgs) poles play a very important role in the calculation of the annihilation cross section; 
and the next-to-lightest neutral exotic L2 can be quite close in mass to L^, in which case 
L'j'L^ co-annihilation as well as L2L2 annihilation have to be considered as well. 

We have followed ref.[^ in our computation of the relic density. For completeness we 
list, but do not derive, the relevant expressions here. The calculation proceeds in two steps. 
First one has to determine the freeze-out temperature Tf, where the relic particles fall out 
of thermal equilibrium. It is given by 

0.0385fefrMp|mi|(crefrw)(x/) 
^/ = 1^ 7= ' (30) 
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where we have introduced Xf = \mi\/Tf, mi being the mass of the hghtest neutral exotic. 
Mp = 1.22 • 10^^ GeV is the Planck mass, and g^, is the number of degrees of freedom that 
are in thermal equilibrium at temperature Tj; we have for simplicity used a fixed ^/g^ = 9, 
which introduces a negligible error in our calculation. Eg. (|30|) has to be solved iteratively, 
since the r.h.s. depends on Xf both explicitly and via the thermal average (acsv), defined as 



{acsv){x) 



X 



3/2 



2v^ 
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(31) 



where v is the relative velocity of the two annihilating particles in their center-of-mass frame. 
In the presence of co-annihilation, the effective annihilation cross section (Teff is given by 



au + 2ai2(l + A)3/2e--^ + ^22(1 + Afe' 



-2xA 



where aij = a{L^L^ anything), and we have introduced 

A 
ges 



- 1; 

mi 

\ + + A)3/2e--^ 



(32) 

(33a) 
(33b) 



Qes is the effective number of degrees of freedom of CDM particles {g = 2 for a single 
Majorana fermion). Numerically, ~ 20 or so for models leading to an acceptable relic 
density. 

Once Xf has been determined, the relic density is given by 



1.07-109 GeV"^ 



J{xf)^Mp 
where the annihilation integral J is defined as 



{acsv){x] 



dx. 



x^ 



(34) 



(35) 



As usual, we have expressed the relic density Q in units of the critical (closure) density, so 
that r2 = 1 corresponds to a flat universe as predicted by inflationary models |jll|. Finally, h 
is the present Hubble parameter in units of 100 km/ (s-Mpc); it lies in the range 0.4 < h < 1. 
The constraint that the universe be older than 10 billion years implies Qh"^ < 1; the true 
upper bound is almost certainly tighter than this, but we want to be conservative since our 
calculation will only be precise on the 10% level. 

In our calculation of the annihilation cross section of eq.(|3^ we have only included 
annihilation into // pairs, where / is a fermion contained in the Standard Model. This 
under-estimates the total annihilation cross section, and hence over-estimates the relic den- 
sity, if is heavy enough to annihilate into pairs of gauge or Higgs bosons, or mixed 
gauge-Higgs final states. However, in this case annihilation into // final states is by itself 
usually sufficient to give an acceptable relic density. The reason is that a large \mi\ also 
implies fairly large SU{2) doublet components of L?, and hence sizable couplings to gauge 
and Higgs bosons. 
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The process L^L^j — * // can proceed via the exchange of a neutral (scalar or pseudo- 
scalar) Higgs boson, or a Z boson, in the s— channel. As discussed in Sec. 3, there might 
also be terms in the superpotential that couple exotic letpons to (s)fermions present in the 
MSSM. However, the resulting contributions to the annihilation cross section are strongly 
suppressed. We already saw that bounds on rare processes severely constrain these cou- 
plings []TU[. Moreover, they would contribute to annihilation only through the exchange of 



a sfermion in the t— or m— channel. The experimental lower bound on squark masses is 



already quite high pS|; furthermore, as mentioned earlier, one expects sfermion masses to 
be of order Mz' in this model ^j, which leads to a strong suppression of sfermion ex- 
change contributions. Finally, these couplings involve the SU{2) doublet components of 
which are usually quite small for i = 1,2; the corresponding contribution to the annihila- 
tion matrix element therefore involves two small mixing factors. This is also true for the 
Z— exchange contribution, but the HL^L^ coupling (|2^ ) contains only one small mixing an- 
gle if the Higgs boson is mostly an SU{2) doublet.| We therefore neglect possible t— channel 
exchange diagrams. 

Including the s— channel exchange of Z and Higgs bosons, the relevant matrix element 
can be written as 

A {L^{h)L%h) f{pi)fip2)) = Af + Af + 4"), (36) 

with 

^ = '^«y'c</).(b),.,.(fa) ^ _ ; „„^r/ (P.)7...(P.); (37b) 

3 A^)Ai) 

- ^^i^M^^M^ g (370 

The ZL^L^ couplings tij are given in eq.(0), while the AL'^L^ couplings t[f^ and the h^L^L^ 
couplings t[jl can be read off eq.(^). The Aff couplings c*^-^^ and hlff couplings are 
listed in eqs.( |194|19b|) . Finally, the Zff couplings a/ and bf are as usual given by a/ = 



— ^li+q^ sin^ 9w, bf = ^/l , with // and q-^ the weak isospin and electric charge of fermion /, 
respectively. In eq. ( P7a| ) we have introduced the total momentum = {ki+k2)i_i = {pi+P2)ii, 
and s = Pf'P^. 

The A and Z exchange diagrams interfere with each other, but not with the scalar Higgs 
exchange contribution. The annihilation cross section can most readily be computed using 
standard trace techniques. The result can be written as 



■l-Recall that this couphng resuhs from the superpotential which couples one singlet to two doublets. 
An SU{2) doublet Higgs boson therefore couples to L^L'j via one singlet and one doublet component. Recall 
also that in this model all light Higgs bosons must be dominantly doublets, since x ':$> v,v in eq.@, so that 
the singlet Higgs boson essentially decouples from physics at scale Mz- 
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where I3f = Jl — Arnj/s, Nf = 1 (3) for leptons (quarks), and the scaled annihilation cross 



sections aij are given by: 



a. 



(ZZ) 



cr, 



(ZA) 



a, 



(AA) 



cr, 



(hh) 



cos9w 



{U,f {a) + b]) 



u jmirrij 



[s-MlY + MlTl 



2m) 
'Ml 



1 - 



2M| ; (s - M|)2 + M|r| 



+ 



V cos6' 



w 



tijt'lf^c^^^bfmf{mi + rrij) 



s — m, 



m'j{mi + rrij)'' 
Ml 



s — (mi — m, 



Ml 



COS^VK 



cos6'vi/ 



) {s-m\f + m\V\' 

3 



s — (mj + m,; 



s(3 



E 

A;=l 



Ah)M) 



mlo +imhoThO 



(39a) 
(39b) 
(39c) 

(39d) 



Here, the initial 3-momentum fc^ is again given by eg. (P^ with M| replaced by s. As before, 
the masses m^ and can have either sign. 

A few comments are in order. When writing the thermal average (^Tj), we have used 
non-relativistic kinematics; for consistency we therefore also have to use a non-relativistic 
expression for s in eqs. (|39|) . 



(|mj| + |mj|)^ + \mimj\v'^. 



(40) 



this might seem dangerous, since in the presence of narrow poles the integral in eg. (pl|) 
can receive sizable contributions from f ~ 1. However, we have checked that using fully 
relativistic kinematics everywhere does not change the result significantly; on the other 
hand, combining eg.(|3TD with a relativistic expression for s can under-estimate the relic 
density by a factor of 2 or 3. 

In principle we now could compute Qh'^ numerically, by inserting egs. (|38D and (|39D into 
egs.(|5nD-(P5D. Note, however, that inserting eg.(PTD into eg.(P5D leads to a double integration. 
Since we want to test several million combinations of model parameters, in order to make 
sure that we covered all relevant regions of parameter space, a direct numerical integration 
is not practical. We used the following approximate method instead. 

As well known |Tl|, |2^, the integral in egs.(^ID and (^) can be computed guite reliably 
from a simple analytical expression if the annihilation cross section aes does not depend too 
sensitively on v. In this case one can use the Taylor expansion 



(41) 



which gives {aesv){x) = + 6bes/x, and J{xf) = aes/xj + 3bcs/x'j. Since Xf ~ 20, annihi- 
lation from an s— wave initial state, which contributes to aes, reduces the relic density more 
efficiently than annihilation from a p— wave does, which only contributes to bes- However, in 
our case this expansion can be used with some reliability only for those contributions that 
do not have an s— channel pole. Specifically, this includes the last term in the sguared Z 
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exchange contribution ( |39a| ), which is due to the longitudinal polarization state of the Z; the 
Z — A interference term (|39b| ) ; and the H\ — iJf* interference terms in eg. ( P9d|) . All other 



contributions do in general show strong variation with f ^, and have to be treated separately. 

Since we use the expansion ( ^ ) for all interference terms, we now only need to compute 
the thermal average over a Breit-Wigner propagator, multiplied with a power of v^: 

/.. . /; C-)^^) ^ (42) 

^V^ "'o \s[v) - mf,] + m^pT^p 

where s{v) is given by eg. (|40D , C is a constant, and P = A oi h\. Setting for simplicity 
s = (I mil + in the numerator (but not the denominator) of eg. (|42|) , /„ can be written 

as 



4 = C^ / dve-^^'/' ^ , , (43) 

2V7r V mirrij J Jo (t;2 - v^) +7 

with 

. = ^^-(KI + KI)'; (44a) 
\mimj\ 

7=f^)'. (44b) 

Note that f q < implies s > m?p for all f > 0, so that the pole is never accessible. On 
the other hand, if fg > 0, the contribution from f ~ t>o to the integral in eg. (|43|) scales like 
g-^''o/'^(t)2)i+"^-o-5^ which can be substantial even for wq ~ 1 if 7 -C 1, i.e. if the pole is very 
narrow. 

The integral in eg. (^3]) still seems to depend on three parameters (the inverse tempera- 
ture, and the position and width of the pole). Further progress can be made by using the 
substitution v' = ^/xv■. 

^2~n / \,yy^ I _i_ . I \ POO .„ { /j/2^1+?^ 





f - vlf + 7 

^Cx'-[ '-' • ) In{vli), (45) 

with -Uq = xfQ and 7 = x^7. We have computed Jq and Ii numerically for 200 values of Wq 
between —100 and +225, and 50 values of 7 between 6.3 ■ 10~^ and 3.5 ■ 10^.| Note that these 
10,000 values of Jq and Ii need to be computed only once; afterwards the thermal average 
over the annihilation cross section can be computed without numerical integration: If Vq 
does not lie in this range, the pole is so distant that the expansion (^) can be used for it; 
for values of Vq and 7 in the specified range, the are estimated by interpolation. 

Specifically, for i = j {L\L\ or L^L^ annihilation), only the sguared A exchange term 
( P9c| ) gets contributions oc Jq when inserted in eg. (|3lD ; the resonant contributions to eg.( P9a|) , 

^Recall that we need to compute /„ only for x > 20. Moreover, Tp/nip > 10~'^ even for the hght Higgs 
bosons of the model. 
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as well as eq.( p9dD , are proportional to v"^ for = m^, i.e. only contribute via Ji. The 
reason is that on-shell Z and scalar Higgs bosons can only be produced from two identical 
Majorana fermions if they are in a p— wave state, while on-shell pseudo-scalar bosons can 
be produced from an s— wave state. However, squared Z and scalar Higgs exchange do 
contribute to /q terms if 7^ m^. In the important special case where L\ and form a 
Dirac fermion, one has mi = —m2] in this case squared A exchange only starts at order f ^ 
(/i terms only), while squared scalar Higgs and Z exchange start at order Finally, for all 
contributions that start at order (Jq terms), we have expanded the numerator to order f^, 
i.e. added a (properly normalized) Ji term. We checked numerically that our combination 
of eq. (|4l|) for non-resonant (interference) terms, and eq. (^) with interpolation to the actual 
values of Vq and 7, reproduces the exact thermal average of the annihilation cross section to 
an accuracy of about 10% or better; this is quite sufficient for us. In contrast, simply using 
the expansion ( ^If ) for the entire cross section can both over- and under-estimate the true 
thermal average by a large factor; this had been observed previously in the similar case 
of the MSSM with scalar Higgs mass close to twice the LSP mass. Finally, the annihilation 
integral (|35|) usually converges rather quickly; we have therefore computed it numerically, 
still using the method outlined above to determine {ae{{v){x) in the integrand. 



5) Results 

We are now in a position to describe our numerical studies of the model, using the expressions 
given in the previous three sections. Our basic approach is to randomly sample the parameter 
space, and to count the number of solutions that satisfy all the constraints we impose, 
including those that can be derived from searches at future (linear) e~^e~ colliders. The 
ultimate goal is to devise a set of constraints such that there are no acceptable solutions left. 
This means that either a signal characteristic for the model has been found, i.e. one of the 
hypothetical future searches is successful, or the model is completely excluded, independent 
of the values of the free parameters. 

We chose this Monte Carlo approach since the number of free parameters is too large 
to allow for a systematic scan of the entire parameter space. The mass matrix (^) for the 
neutral exotic leptons already contains ten free parameters. For each scan of parameter 
space, we impose a lower bound on the charged leptons masses mj^±, since these particles 

i 

should be trivial to discover at an e^e~ collider unless their production is kinematically 
suppressed. We also fix the ratio tan/5 = v/v for a given run; this determines v and v, since 
+ = 2M]v/g is known. The Yukawa couplings appearing in eq.(^) are then chosen 
randomly in the interval \Xijk\ < Amax- We take Amax = 0.85, which is approximately the 
upper bound on any one coupling from the requirement that there should be no Landau 
pole at scales below the GUT scale. This is a conservative approach, since this bound is 
significantly stronger if several Yukawa couplings are sizable |TD|, which is required if m^o is 
to be close to its upper bound. 

Having specified the exotic lepton sector, we check whether the present LEP constraints 
( p5D and ( p6D are satisfied. If so, we check whether the total "visible" exotic cross-section, 
determined from eq . (^7]) , is sufficiently small so that a given e~^e~ collider will not detect 
pair production of exotic leptons; this constraint obviously depends on the collider we are 
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considering, as described below. Recall that we consider a neutral exotic lepton to be 
"visible" only if it is at least 5 GeV heavier than L^. 

This exhausts the constraints from searches for exotic leptons at present and future e~^e~ 
colliders. If our choice of parameters is still viable, we next randomly pick a Higgs sector, 
subject to the constraints (|28|) . Since present bounds on the mass of the Z' boson are already 
quite high [Q, we always find that there is one singlet-like neutral Higgs boson, which plays 
no role in any of the process we are considering; it is therefore sufficient to simply fix the SM 
singlet vev x to some large value, say 2 TeV. Since our runs are for fixed tan/5, we only need 
to chose the values of two additional parameters in order to completely specify the Higgs 
sector of the model, see eq.(|). We chose these to be the coupling A333, which also has to lie 
in the interval IA333I < 0.85, and the mass tua of the physical pseudo-scalar Higgs boson; 
this then fixes the parameter ^333 via eq.(|l^). Note that for each set of leptonic parameters 
that satisfy constraints from e'^e~ colliders, we keep chosing pairs (A333, m^) at random until 
we have found an acceptable Higgs sector. 

This is necessary because Higgs exchange contributions can play an important role in the 
calculation of the L° relic density. Once the exotic lepton sector and the Higgs sector are 
specified, all quantities appearing in the annihilation cross sections (^) are fixed, and we 
can compute the relic density as described in Sec. 4, and check whether it is acceptable. 

Figs. la,b show that present constraints [the bounds (^Sf ), (^) and (^8|) , together with 
Qh"^ < 1] are still quite far from testing the model decisively. For these figures we have 
allowed charged lepton masses as low as 40 GeV when sampling the parameter space; we 
have however required |?7i^±| > 45 GeV for all acceptable solutions. This explains why the 

i 

survival fraction after imposing LEPl constraints only is less than unity in some bins with 
m^o > Mz/2. 

Note that LEPl searches only impose weak constraints for very small values of rrij^o. 
The reason is that most parameter sets with small rrij^o have large masses for the charged 
exotic leptons. In this case four of the six neutral exotics are also heavy, while the remaining 
two states are mostly SM singlets, i.e. couple only weakly to the Z boson. However, the 
relic density constraint is most effective precisely in this situation, since it leads to small 
annihilation cross sections: The couplings of singlet-like exotics to the light Higgs bosons 
are also weak, and their small masses suppress the cross sections even further; far below 
the pole, s— channel exchange contributions scale like m?o/mp, where P is the exchanged 

i 

particle. We therefore find no cosmologically acceptable scenario with m^o < ~ 5 GeV; 
this is not surprising since the bb final state, if accessible, contributes most to Higgs exchange 
diagrams. 

Figs. 1 also reveal a technical problem: even though each figure is based on 10^ sets of 
leptonic parameters ("hundred thousand models"), the region of large rrij^o is only sparsely 
populated. This is perhaps not surprising, since all ten parameters appearing in the mass 
matrix (^) must be chosen within a narrow range if m^^o is to come out close to its upper 
bound. At the same time, combinations of parameters leading to large m^o are most difficult 
to exclude, i.e. most easily evade all bounds. For one thing, large m^o implies large SU (2) 
doublet components of L^, and hence large annihilation cross sections and a small relic 
density. Further, constraints from collider searches can most easily be evaded if is close 
in mass to L°, since then also becomes effectively invisible. In such a situation L^L2 
co-annihilation can also reduce the relic density even further. Since our goal is to test the 
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model decisively, it would be advantageous if in our sampling of parameter space we could 
give preference to regions that are most difficult to exclude. 

Fortunately a large m^o is correlated with a small — mass difference. In ref.[|lO 



it has been shown that choices of parameters that maximize rrii^o always lead to a situation 
where the six eigenvalues of the mass matrix come in three pairs that only differ by 
a sign; in such a situation the six Majorana states can also be described by three neutral 
Dirac fermions. In particular, the Majorana state will now be completely invisible, being 
degenerate in mass with and hence (by assumption) stable; this obviously also maximizes 
L1L2 co-annihilation. From now on we will therefore only show results for the subset of 
parameter space where the six neutral exotic leptons do indeed form three Dirac fermions. 
This can be enforced by chosing 

Al31 = Ai32, A231 = A232, A311 = — A312, A321 = — A322- (46) 

For given tan/?, this reduces the number of parameters in the exotic lepton sector from ten to 
six; this reduced parameter space is obviously much easier to sample exhaustively. We have 
checked that our runs with eqs.(^) imposed do always find significantly more acceptable 
solutions than scans of the entire parameter space with equal statistics. 

Figs. 2a,b show that, at least if the mass matrix (^ has Dirac structure, neutral lepton 
searches at LEP2 will not lead to significant new constraints, unless a light charged exotic is 
found. We have assumed here that LEP2 will reach a center-of-mass energy ^/s of 190 GeV, 
and that the production of neutral leptons is detectable if the cross section, summed over 
all "visible" modes, exceeds 20 fb, which corresponds to 10 events per experiment for the 
foreseen integrated luminosity of 500 pb^^. Comparison with Figs. 1 shows that the fraction 
of parameter space excluded by present LEPl constraints has become smaller, even for light 
This is mostly because L^L2 and L2L2 final states are now invisible. Note also that 
we now find some cosmologically acceptable solutions with tan/3 = 1.2 and m^o as small as 
10 GeV; this indicates that co-annihilation can indeed reduce the relic density significantly. 
Further, even though the absolute upper bound on m^o decreases with increasing mass of the 
charged exotic leptons [|l^. Figs. 2a,b extend to larger values of m^o than Figs la,b do; clearly 



the restrictions (^) have made it much more likely to produce scenarios with large m^o. 
Finally, comparison of Figs. 2a and 2b shows that LEP constraints exclude a larger fraction 
of parameter space for large tan/?; this can already be seen from Figs. la,b. We showed 
in Sec. 2 that increasing tan/? decreases the upper bound on m^o, since detM.LO oc (vvx)'^; 
however, the size of the SU{2) doublet components of the light exotics, and hence their 
couplings to the Z boson, remains more or less the same. Reducing the masses of the exotics 
while leaving their couplings essentially unchanged obviously increases the partial width for 
Z decays into exotic leptons. 

Figs. 2 clearly show that LEP2 will not be able to test the model decisively; there are 
many choices of parameters that lead to an acceptable cosmology, but no "new physics" 
signal at this collider. One will therefore need (linear) colliders operating at higher energies 
in order to probe the entire parameter space. Such colliders are often assumed to be built in 
three stages, where the energy is increased from about 0.5 TeV to 1.0 TeV and, eventually, 
1.5 TeV or even higher. We generically call these three stages NLCl, NLC2 and NLC3. 

Figs. 3a,b show the situation if a 500 GeV e~^e~ collider fails to discover pair production 
of exotic leptons, which we have interpreted as meaning that the total cross section into 



18 



visible exotic final states is less than 0.5 fb; the integrated luminosity of such a collider is 
usually assumed to be several tens of fb~^ per year. For these runs we have required the 
masses of the charged exotic leptons to exceed 240 GeV. This reduces both the upper bound 
on rrij^o and the maximal SU{2) doublet component of significantly. As a result, for 
tan/3 = 1.2 present LEPl searches do not constrain the model any further. For tan/3 = 5, 
these constraints still do exclude some combinations of parameters, but they are clearly 
much less restrictive here than for lower masses of the charged exotics, see Fig. 2b. On 
the other hand, unlike at LEP2, searches for neutral exotics at NLCl can probe sizable 
regions of parameter space even if no charged exotic leptons are found; this can be seen from 
the large differences between the dotted and dashed histograms in Figs. 3. Finally, the relic 
density constraint again excludes parameter choices giving m^^o < rrib, and impose significant 
constraints as long as rrij^o < 15 to 20 GeV. However, a substantial number of parameter 
sets still satisfies all constraints (solid histograms). 

This remains true even if a 1 TeV collider fails to discover pair production of exotic 
leptons. Indeed, Figs. 4a,b show that a small region of parameter space survives even if 
searches for pairs of exotic leptons at a 1.5 TeV collider remain unsucessful; here we have 
assumed that the integrated luminosity scales like the square of the beam energy, so that 
a signal exceeding 0.05 fb would be detectable. We evidently need to find some additional 
constraint (s) if we want to test the model decisively. 

There are significant differences between Figs. 4a,b and 3a,b, which might offer a clue as 
to what these additional constraints might be. To begin with, the minimal allowed value of 
m^o for small tan/? has increased from about 6 GeV (Fig. 3a) to about 23 GeV (Fig. 4a). 
Unsuccessful earches at a 1.5 TeV collider force the SU{2) doublet components of to be 
much smaller than searches at an 0.5 TeV collider do, which leads to considerably reduced 
couplings of the light exotic leptons to gauge and Higgs bosons. This has to be compensated 
by an increase of m^o in order to keep the relic density acceptably small; recall that for 
light L°, the annihilation cross sections scale like m^o. Even more importantly, if a 1.5 TeV 
collider fails to find evidence for the pair production of exotic leptons, the upper bound on 
m^o is reduced to about 30 (11) GeV for tan/3 = 1.2 (5), as compared to 84 (27) GeV after 
unsuccessful searches at NLCl. 

Recall that the masses of the two lightest exotic Majorana states, or the lightest exotic 
Dirac state if eqs. (^6]) hold, come from vevs that break SU{2). This indicates that there 
is a correlation between m^o and the size of the couplings between light exotics and SU (2) 
doublet Higgs bosons, although the relation is not as simple as that between the masses and 
Yukawa couplings of the quarks and leptons of the SM. Further, Higgs searches at LEPl 
imply that the mass of a scalar Higgs boson with unsuppressed ZZH coupling must exceed 
60 GeV, see eg. (|28a|) . Such a Higgs boson can therefore decay into pairs of light neutral 
exotic leptons for all allowed combinations of parameters shown in Figs. 4; this could lead 
to a large invisible branching ratio of this Higgs boson. Indeed, we find that all surviving 
scenarios shown in Fig. 4a have one rather light neutral scalar Higgs boson, with essentially 
unsuppressed coupling to two Z bosons, and with invisible branching ratio exceeding 80%. 
Moreover, it will be produced copiously at any e+e~ collider with ^/s > 300 GeV, since 
its mass cannot exceed 150 GeV or so, see eq. (|T0|) . Note that such a Higgs boson is easily 
detectable even if it decays invisibly [^, since it would be produced in association with a Z 
boson, whose decay products would be sufficient to reconstruct rriH- Finally, a large invisible 
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branching ratio for a hght scalar Higgs boson cannot be accommodated in the SM; in the 
MSSM it would imply the existence of a chargino light enough to be discovered at LEP2 



2J]. It therefore constitutes a signal for the kind of model we are considering. 

The situation at large tan/?, Fig. 4b, is somewhat more complicated. Notice that, unlike 
in Fig. 4a, now only a small fraction of all parameter sets that give m^o close to its upper 
bound, and do not lead to a detectable signal for exotic lepton pair production, survives the 
relic density constraint. The reason is that for small tan/5 and with m^o close to its upper 
bound, Z exchange by itself is usually sufficient to give an acceptable relic density. On the 
other hand, increasing tan/? to 5 reduces the maximal m^o so much that Z exchange, the 
cross section for which scales like m^o/M|, is no longer sufficient; light Higgs bosons have 
to be present to enhance the annihilation cross section. In particular, a light pseudo-scalar 
Higgs boson A allows L^L^ and L2L2 to annihilate from an s— wave initial state, which gives 
a larger thermal average than annihilation from a p— wave initial state does, as discussed 
below eq.(^) in Sec. 4.Q Therefore even if rriio is close to its upper bound, the relic density 
will only be acceptable if rriA is chosen to be fairly small, which is true only for some fraction 
of Higgs parameter space. We find that all surviving parameter sets in Fig. 4b have niA < 110 
GeV, and also have two neutral scalar Higgs bosons with masses below 125 GeV. At least 
one of these three light neutral Higgs bosons has invisible branching ratio exceeding 50%. All 
three of these Higgs bosons will be produced copiously, either in association with a Z boson 
or as Ah^ pairs. A large invisible branching ratio for any of these Higgs particles is therefore 
again a distinctive signature. Note that now the invisible branching ratio is always less than 
75%, so that in at least 50% of all Ah^ pairs at least one of the two Higgs bosons will decay 
into a visible final state, mostly bb pairs. This not only guarantees the detectability of this 
final state, but also ensures that the invisible branching ratios can be measured with some 
accuracy, e.g. by comparing the rate for events with one invisible Higgs boson (single-sided 
events) with that for events where both Higgs particles leave a detectable final state. Finally, 
we note that in the surviving cases the charged Higgs boson mass (0) also is below 125 GeV, 
and usually below 100 GeV; such a light charged Higgs boson is already almost excluded in 
the MSSM, and might thus give a second indication for physics beyond the MSSM. 

Based on the results of Figs. 4 we therefore conclude that in this model either a signal for 
the pair production of exotic leptons will be observed, at the latest at an e"'"e~ collider with 
^/s = 1.5 TeV; or at least one neutral Higgs boson with mass below 150 GeV and invisible 
branching ratio exceeding 50% will be found at any e~^e~ collider with s/s > 300 GeV. If 
neither of these two signals is detected, the model can be completely excluded. 

Since this conclusion is based on a large but finite Monte Carlo sampling of parameter 
space (Figs. 4a,b contain 5 ■ 10^ sets of leptonic parameters each), we have checked that it 
also holds in a "worst case" scenario. To this end we have chosen the parameters entering 
the mass matrix (^ such that rrij^o is maximized |TD]: 



Al31 — Ai32 — A311 — — A312 — A23I — — A321 — A232 — A322 — Amax! (47a) 

"^Lf = -"^L±- (47b) 



*L]'I/2 co-annihilation via scalar Higgs exchange could come from an s— wave initial state; however, if 
the neutral lepton mass matrix has Dirac structure, the off-diagonal L^L'^h'^ couplings vanish identically, as 
does the LfL^A coupling. 
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Note that this ansatz is consistent with eqs.(^). In Fig. 5 we show the resulting maximal 
value of m^o as a function of tan/?, subject to the constraints > m^^i.min and Imiol + 
\mio\ > 2m2,±,min- The latter bound approximates the constraint that can be derived from 
an unsuccessful search for associate LlL'^ production; note that the relevant coupling is 
always quite large for the ansatz (^TD-H For our choice Amax = 0.85, the maximal allowed 
m^o for mL±^rnin = 700 GeV is indeed around 30 GeV, in agreement with results shown in 
Fig. 4a; for tan/3 = 5, m^o falls to about 11 GeV, in agreement with Fig. 4b. We should 
mention here that the choice Amax = 0.85 in eq.( [47a| ) is very conservative. In ref.|]lO| it has 
been pointed out that Amax > 0.7 implies the existence of a Landau pole at a scale below 
lO^'' GeV, the smallest energy scale where one might expect the gauge group (and hence the 
relevant renormalization group equations) to change in this class of models. Note that for 
large m^i mim the upper bound on m^o scales like the square of Amax; reducing Amax to 0.7 
therefore means that m^^o < 30 GeV already for mi±^min = 480 GeV, in which case a 1 TeV 
collider would be sufficient to test the model decisively in the manner described above.|^ 

We also checked whether it is possible to chose parameters such that the SM-like Higgs 
boson only has very weak couplings to light neutral exotics. This is indeed possible, but 
only for sizable tan/?, and only if the pseudoscalar Higgs boson is quite light. In this case the 
model contains two light SU{2) doublet neutral scalar Higgs bosons, as mentioned earlier; 
the "SM-like Higgs" is defined to be the one with the larger coupling to two Z bosons. We 
saw in Fig. 4b that for large tan/3, the relic density constraint can only be satisfied if rriA is 
rather small; the conditions for a scenario with small invisible width of the SM-like Higgs are 
therefore satisfied. However, Fig. 6 shows that either the pseudoscalar Higgs or the lightest 
neutral scalar can always decay into an L^L^ pair if searches at a 1.5 TeV collider do not find 
a signal for the pair production of exotic leptons. In this figure we show the minimal allowed 
values of and of m^^o, as well as the minimum of the sum of these two masses, as a function 
of tan/?, where we have only used the present LEPl constraints (pH]). For tan/3 ^ 1.5, 
could be as low as 34 GeV, so that A — >■ L^L'l decays would be kinematically forbidden over a 
sizable region of parameter space. However, in this case the lightest scalar Higgs boson must 
be heavier than 50 GeV, since rriA + > 84 GeV. Note that in this scenario, /i° is usually 
not the SM-like Higgs, and does have sizable couplings to light exotic leptons. Moreover, we 
find that the invisible branching ratio of the SM-like scalar Higgs can only be reduced to a 
value below 10% if tan/5 > 3, in which case unsuccessful searches for lepton pair production 
at NLC3 imply m^o < 20 GeV, see Fig. 5. In this case both A and hi have large invisible 
branching ratios; remember that a sizable contribution from A exchange to annihilation, 
and hence a substantial AL^L^ coupling, is required to satisfy the relic density constraint for 
tan/5 > 2. Finally, we found no acceptable solutions where the invisible branching ratio of a 
light Higgs boson can be diluted by decays into pairs of even lighter Higgs bosons {h^ hlhl 
or AA). We conclude that it is indeed impossible to devise a Higgs sector such that 

no light Higgs boson has large invisible branching ratio if NLC3 fails to find a signal for the 
pair production of exotic leptons. 

Finally, we investigated the question whether Higgs searches at LEP2 will allow us to 

^Recall that the production of two hght neutral exotic leptons leads to an invisible final state. 
^The larger value of A,nax we are using here also explains why Fig. 5 allows larger values of mj^o than 
ref.0 does. 
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sharpen our predictions. This does not seem to be the case, at least as far as searches at a 1.5 
TeV coUider are concerned. As mentioned earher, for small tan/3, Z exchange by itself can 
be sufficient to give an acceptable relic density, so increasing the experimental lower bounds 
on Higgs masses will have little effect on Fig. 4a. For tan/? > 3 one needs vtia < HO GeV in 
order to achieve < 1. This upper bound is almost independent of tan/5 once it exceeds 
3 or so; the decrease of m^o caused by increasing tan/5 is balanced almost perfectly by the 
increase of the Abb and At^t~ couplings, leading to a constant A exchange contribution, as 
long as m^o > m;,. However, this bound still allows values of just beyond the reach of 
LEP2. 

On the other hand, if LEP2 can increase the lower bound on the SM Higgs to 90 GeV or 
more, a 1 TeV e^e~ collider would be sufficient to test the model decisively even with our 
choice Amax = 0.85, since the absence of a signal for exotic lepton pair production at such a 
collider implies m^o < 45 GeV, see Fig. 5; one would then again be in the situation where 
at least one copiously produced Higgs boson has to have a large invisible branching ratio. 
Unfortunately in the absence of a positive signal for lepton pair production, a 0.5 TeV collider 
would never be sufficient to test the model uniquely, since, as shown in Fig. 3a, if would 
only allow to establish the bound mjj^^ < 85 GeV, well above half the maximal allowed mass 
of the SM-like Higgs boson. One can then always find cosmologically acceptable scenarios 
where only one Higgs boson is experimentally accessible, with small or vanishing invisible 
branching ratio; such a model would be indistinguishable from the MSSM, or even the SM, 
as far as the NLCl is concerned. 

6) Summary and Conclusions 

In this paper we have shown that searches for exotic leptons and Higgs bosons at future 
e+e~ coUiders, when combined with cosmological constraints, can test a large class of £'(6) 
models decisively. The basic idea is that, due to the structure of the mass matrix for the 6 
neutral exotic leptons predicted by the model, the failure to detect the production of heavy 
(charged or neutral) exotic leptons at such colliders will allow to derive increasingly stronger 
upper bounds on the masses of the fight neutral exotics. Eventually these leptons will have 
to be lighter than half the mass of the light neutral Higgs bosons of the model. The only way 
to suppress the couplings of these light exotic leptons to all light Higgs bosons is to make 
the leptons almost pure SU (2) singlets. However, in this case exotic leptons produced in the 
very early universe would have small annihilation cross sections, so that their present relic 
density would be unacceptably large. All cosmologically acceptable parameter sets therefore 
predict that either the pair production of exotic leptons is visible at a 1.5 TeV e'^e~ collider, 
or at least one light neutral Higgs boson will have an invisible branching ratio exceeding 
50%; if LEP2 fails to discover a Higgs boson, a 1 TeV e^e~ collider would be sufficient. 

In order to arrive at this conclusion, we had to make two crucial assumptions. First, no 
5*0(10) singlet scalar N is allowed to have a vev at an "intermediate scale" of O(10^°) GeV 
or more. This assumption is not unreasonable, since such a large vev would allow to push 
the masses of almost all new particles predicted by £"(6) up to this high scale, including 
the new gauge bosons, the exotic quarks and leptons, and their superpartners. Apart from 
the possible presence of light right-handed (s)neutrinos, the model would then look like the 
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MSSM at scales below this vev. Note that we do allow 5*0(10) nonsinglets to have such a 
large vev. Secondly, we have to assume that i?— parity is conserved. Otherwise the lightest 
exotic lepton, which has odd i?— parity, could decay even if it is lighter than the lightest 
neutralino, in which case the relic density constraint would be satisfied trivially. However, if 
-R— parity is broken, this rather complicated model would not be able to accommodate cold 
dark matter; recent attempts to understand structure formation in the universe strongly 
favor scenarios with a substantial amount of cold dark matter ||26[| . 

We also made a few additional assumptions in order to simplify our calculation. However, 
since we have been very conservative in our interpretation of the cosmo logical constraint, 
and of the upper bound on the Yukawa couplings of the exotic leptons that follows from the 
requirement that no coupling should have a Landau pole below the intermediate scale of 10^° 
GeV, we believe that our final result holds even if these additional, technical assumptions 
are relaxed. Specifically, we have assumed that the exotic leptons do not mix with charginos 
and neutralinos. This assumption is technically natural. Moreover, in ref. it has been 
argued that, while such mixing might increase the upper bound on 771^0 for fixed Yukawa 
couplings, it would also necessitate the existence of additional terms in the superpotential; 
these would lower the upper bound on the relevant Yukawa couplings from the absence of 
Landau poles. The total change of the upper bound on m^o is therefore quite modest. 

Further, when computing the annihilation cross sections, we have assumed that there 
are no terms in the superpotential which couple exotic leptons to ordinary quark and lepton 
superfields. However, we argued in Sec. 4 that such couplings would in any case contribute 
negligibly to the total annihilation cross section. Finally, we have taken the Z' boson(s) of 
the model to be too heavy to contribute either to the pair production or to the annihilation of 
exotic leptons. Present lower bounds on Mz' already ensure that Z' exchange contributions 
to the annihilation cross sections are negligible. They could still affect the production of 
exotic leptons at TeV-scale e+e~ colliders, but in this case the existence of a Z' boson could 
also be inferred from studies of ordinary quark and lepton pair production, which would 
again give a good signal for the class of models we are studying. Besides, the constraints on 
the exotic lepton sector that could be derived from new particle searches at future colliders 
are likely to be more, not less, severe if there is a significant contribution from Z' exchange, 
since even the SU{2) singlet fields N couple to the Z' boson with full gauge strength. 

There is one more assumption that we have not mentioned so far: when using the for- 
malism described in Sec. 4 to estimate the relic density, we have assumed that the exotic 
leptons are non-relativist ic ("cold") when they drop out of thermal equilibrium. This is 
true if their masses exceed a few hundred MeV or so, and the relic density constraint will 
remain valid down to much smaller masses, in the keV range, but our calculation clearly 
breaks down if the light exotics are (nearly) massless. This could, for example, be achieved 
by setting all couplings in the superpotential (P to zero, except for A113 and A223 which 
are needed to give masses to the charged exotic leptons. The two lightest neutral exotics 
would then be massless SU{2) x U{1)y singlets. However, this would increase the density 
of relativistic particles, and hence the expansion rate of the universe, in the epoch when 



light nuclei are formed. The most recent analysis |^ finds that data seem to favor models 
where the effective number of SM neutrinos is smaller, not larger, than three. A model with 
additional massless fermions is therefore strongly disfavored. 

Finally, it can be argued that the upper bound (0) on the mass of the lightest neutral 
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scalar Higgs boson of the model allows a much easier test. However, a very similar bound 
also holds in the MSSM and, indeed, in all SUSY models where the Higgs sector is required 
to remain perturbative up to some high scale [^2|, ITSf. This test would therefore not be 



very specific. While the failure to detect such a Higgs boson would rule out a very large 
class of models, including the E{Q) models we are considering here, discovery of the Higgs 
boson may not be sufficient to distinguish between the present model and the MSSM, or 
even the non-supersymmetric SM. In contrast, the large mass splitting between the light and 
heavy exotic leptons implies that exotic lepton production should not be confused with the 
production of charginos and neutralinos predicted by the MSSM. Large invisible branching 
ratios for light Higgs bosons also clearly indicate the presence of (super) fields beyond those 
contained in the (MS)SM, unless the mass of the light chargino is close to its present lower 



bound. Singlet Majoron models ||33| can also have light Higgs bosons with large invisible 
branching ratios |Q. However, in such models the light bosons have sizable SU(2) singlet 
components, in contrast to the models we are studying here, where the large mass of the Z' 
boson forces all light Higgs bosons to be predominantly SU{2) doublets. A detailed study of 
the production and decay of the light Higgs boson(s) should therefore be able to distinguish 
between the E{6) model and the singlet Majoron model. 

Finally, we would like to emphasize that the decisive test we have devised here relies on 
the versatility of high energy e~^e~ colliders. Clearly the cross sections we are studying are 
too small to give viable signals at hadron colliders. Moreover, the search for the production of 
neutral leptons, in particular the associate production of a light and a heavy neutral exotic, 
plays an important role in our analysis. This would not be feasible at 67 or 77 colliders, 
which may also be unable to discover, let alone study, Higgs bosons with large invisible 
branching ratios. All these colliders will be able to impose constraints on, or - with luck - 
to discover, some of the new particles predicted by E{6) models. However, only searches at 
e"'"e~ colliders, when combined with cosmological considerations, seem capable of excluding 
these models completely. 
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Figure Captions 



Fig. 1 The fraction of parameter sets that survive certain constraints is shown as a function 
of the mass m^o of the hghtest neutral exotic lepton, for tan/5 = 1.2 (a) and 5 (b). The 
fraction is relative to the total number of sets generated in a given bin. For the dotted 
histograms only the current, LEPl constraints (|25|) , (p6|) have been imposed, while the 
solid histograms show the fraction of parameter sets that in addition satisfy Qh'^ < 1. 
Recall that we only include // final states in the calculation of the relic density; this 
under-estimates the annihilation cross sections, i.e. over-estimates the relic density, 
for m^o > rriw- When chosing the parameters of the neutral lepton mass matrix (|^), 
we have required the charged exotic leptons to be heavier than 40 GeV. 

Fig. 2 The fraction of parameter sets that survive present LEPl constraints (dotted), exotic 
lepton searches at LEP2 (dashed), and the relic density bound Qh"^ < 1 (solid). We 
have assumed that LEP2 will be able to detect pair production of neutral exotic leptons 
if the total cross section for final states where at least one leptons is more than 5 GeV 
heavier than exceeds 20 fb. For reasons explained in the text, we have imposed 
the restrictions (|iB| ) on the parameters of the neutral lepton mass matrix (^), which 
implies that the 6 neutral Majorana leptons pair up to form 3 Dirac states. We have 
required the charged exotic leptons to be heavier than 80 GeV here. Results are for 
tan/3 = 1.2 (a) and 5 (b). 

Fig. 3 As in Fig. 2, except that we have required the charged exotic leptons to be heavier 
than 240 GeV, and the dashed histograms show the fraction of parameter sets that 
would not lead to an observable signal for the pair production of exotic leptons at a 
500 GeV e~^e~ collider ("NLCl"). Here a signal is considered to be observable if the 
total cross section for the production of neutral exotic leptons, with one final state 
lepton being at least 5 GeV heavier than L^, exceeds 0.5 fb. 

Fig. 4 As in Fig. 2, except that we have required the charged exotic leptons to be heavier 
than 700 GeV, and the dashed histograms show the fraction of parameter sets that 
would not lead to an observable signal for the pair production of exotic leptons at a 1.5 
TeV e^e~ collider ("NLC3"). Here a signal is considered to be observable if the total 
cross section for the production of neutral exotic leptons, with one final state lepton 
being at least 5 GeV heavier than L^, exceeds 0.05 fb. As explained in the text, all 
surviving parameter sets (solid histograms) predict a large invisible branchig ratio for 
at least one light neutral Higgs boson, which allows to test the model decisively. 

Fig. 5 The maximal value of the mass rrij^o of the lightest neutral exotic lepton, computed 
from the ansatx (|^) with Amax = 0.85, is shown as a function of tan/3. We have 
required m^± > rrii+^^i^, and \mj^o \ + \mio\ > 2m/^+ for different values of rrii+^^i^ 
as indicated. 
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6 The minimal masses of the Higgs bosons of the model, as determined from present 
LEPl constraints p8|), is shown as a function of tan/?. The solid, dashed and dotted 
curves show the smallest allowed mass of the lightest neutral Higgs scalar, of the 
pseudoscalar, and the minimal sum of the masses of the lightest scalar and pseudoscalar 
Higgs bosons, respectively. The minimum of the sum is significantly larger than the sum 
of the minima, indicating that both Higgs masses cannot be minimized simultaneously. 
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